Draft version January 26, 2013 

Preprint typeset using LAT^X style cmulatcapj v. 11/10/09 



O 

(N 



O 

m 

O 

u 

43 

Q-i! 
6 



> 

00 

in 
o 



13 



COSMIC RAY DIFFUSION FRONTS IN THE VIRGO CLUSTER 
William G. Mathews 1 and Fulai Guo 1 

Draft version January 26, 2013 

ABSTRACT 

The pair of large radio lobes in the Virgo cluster, each about 23 kpc in radius, have curiously sharp 
outer edges where the radio-synchrotron continuum flux declines abruptly. However, just adjacent 
to this sharp transition, the radio flux increases. This radio limb-brightening is observed over at 
least half of the perimeter of both lobes. We describe slowly propagating steady state diffusion fronts 
that explain these counterintuitive features. Because of the natural buoyancy of radio lobes, the 
magnetic field is largely tangent to the lobe boundary, an alignment that polarizes the radio emission 
and dramatically reduces the diffusion coefficient of relativistic electrons. As cosmic ray electrons 
diffuse slowly into the cluster gas, the local magnetic field and gas density are reduced as gas flows 
back toward the radio lobe. Radio emission peaks can occur because the synchrotron emissivity 
increases with magnetic field and then decreases with the density of non-thermal electrons. A detailed 
comparison of steady diffusion fronts with quantitative radio observations may reveal information 
about the spatial variation of magnetic fields and the diffusion coefficient of relativistic electrons. On 
larger scales, some reduction of the gas density inside the Virgo lobes due to cosmic ray pressure must 
occur and may be measurable. Such X-ray observations could reveal important information about 
the presence of otherwise unobservable non-thermal components such as relativistic electrons of low 
energy or proton cosmic rays. 

Subject headings: cosmic rays; diffusion; radio continuum; galaxies: clusters; galaxies: clusters: indi- 
vidual Virgo; magnetic fields 



1. INTRODUCTION 

Several nearby galaxy clusters contain extended non- 
thermal radio lobes that appear to be unrelated to cur- 
rently visible X-ray cavities. For example, the iconic 
image of the Virgo cluster at 90 cm (325 MHz) from 
Owen, Eilek & Kassim (2000), shown in Figure 1, fea- 
tures two large opposing quasi-circular radio lobes. (The 
famous 2 kpc radio jet in M87 is over-exposed in the 
central dark region.) In Mathews & Brighenti (2008a) 
we proposed that the Virgo radio lobes were formed by 
relativistic electrons that diffused through the walls of 
X-ray cavities that have long-since buoyantly risen and 
disappeared from view. This conjecture seemed attrac- 
tive because a 25-kpc long radial thermal X-ray filament 
(Young, Wilson & Mundell 2002; Forman, et al. 2005; 
Werner, N. et al. 2010; Million et al. 2010) passes in 
projection right through the center of the southern Virgo 
lobe, suggesting that both were formed in the same X- 
ray cavity event. Moreover, the radio-synchrotron spec- 
tral age of the oldest and most distant electrons in the 
Virgo lobes, t sync m 10 s yrs, matches the dynamical 
age of the post-cavity filament as computed by Math- 
ews & Brighenti (2008a) who suggested that the radio- 
synchrotron electrons currently filling the southern lobe 
may have diffused through the walls of the same filament- 
forming X-ray cavity. Cosmic ray electrons freely dif- 
fusing isotropically from the center of the southern lobe 
of radius r; (, e m 23 kpc suggests a diffusion coefficient 
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filament and southern Virgo lobe is the absence of a cor- 
responding radial thermal filament in the northern lobe. 

In addition, we now suspect that the Virgo radio lobes 
have a more complex evolution resulting from a mul- 
titude of active galactic nucleus (AGN) outbursts that 
continuously resupply cosmic ray electrons to the lobes, 
a point of view also suggested by Owen, Eilek & Kas- 
sim. The 8-shaped Virgo lobes may represent an approx- 
imately bipolar energy ejection pattern from the central 
black hole in M87. AGN jet and cavity activity in Virgo 
(as in Perseus and many other clusters) rarely occurs 
along a fixed bipolar axis as defined for example by an 
unchanging black hole spin, but jet energy may be dis- 
tributed in directions near such an axis. As a conse- 
quence, the Virgo lobes may expand much slower and be 
older than implied by free cosmic ray diffusion. 

With this in mind, we draw attention in Figure 1 to the 
sharp transition at the outer boundary of the Virgo radio 
lobes where the radio flux abruptly decreases, completely 
unlike the smooth gaussian transition expected for free 
diffusion with k ~ 10 30 cm 2 s _1 . Owen, Eilek and Kas- 
sim (2000) emphasize that the sudden decrease in ra- 
dio flux at the lobe boundary is real, not an artifact of 
image-processing software. Another curious widespread 
feature just at this lobe-cluster gas transition are the 
bright synchrotron-emitting rims visible in Figure 1 im- 
mediately adjacent to the sudden sharp drop to zero ra- 
dio flux. According to Owen, Eilek and Kassim, this 
counterintuitive radio limb-brightening is present in over 
half of the perimeter of both lobes. 

In the following discussion we propose that the bright 
radio rims can be understood if the diffusion coefficient 
of synchrotron-emitting electrons decreases sharply near 
the lobe-cluster gas boundaries. 
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2. PLANE PARALLEL STEADY DIFFUSION FRONT 

We assume that cosmic rays (relativistic electrons 
and/or protons) in the radio lobes are coupled to the 
gas by means of a magnetic field frozen into the gas 
through which cosmic rays may diffuse. While the mag- 
netic field allows cosmic ray pressure gradients to act 
on the cluster gas, the small magnetic energy density is 
dynamically negligible. Magnetic fields of strength 0.3- 
10 pG are ubiquitous in cluster gas (Govoni & Feretti 
2004). In general the magnetic field decreases with clus- 
ter radius and decreasing gas density. An approximate 
powerlaw relation B oc n% with 0.4 < s < 0.7 in the 
Coma cluster gas has been established from Faraday de- 
polarization observations of cluster galaxies (Bonafcde 
et al. 2010). However, fields with energy densities 
B 2 /8ir = 6x 10" 13 (B/4/xG) 2 erg cm" 3 cannot signifi- 
cantly influence the dynamics of local cluster gas with 
a much larger thermal energy density 3P/2 = 5 x 
10- n (n e /0.01 cm- 2 )(T/keV) erg cm" 3 , and B 2 /8tt « 
3P/2 probably holds at every radius in the cluster gas. 
In view of this, magnetic fields need not explicitly appear 
in the dynamical equations: 



dp 

dt 
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where e = P/( 7 - 1) = 3P/2 and e c = P c /( 7c - 1) = 
3P C are the energy density-pressure relations for gas and 
cosmic rays respectively and k is the difussivity of cosmic 
rays. Weak magnetic fields evolve passively, advecting 
with the gas. These equations are discussed further in 
Mathews & Brighenti (2008a,b). 

Recognizing the relatively short time for gas to flow 
past the radio lobe boundary, we do not consider the 
energy lost by thermal or non-thermal radiation as gas 
flows through the boundary rims. The time required 
for gas flowing at velocity u w 20 km s _1 to cross 
the bright ratio rims of thickness Ax < 0.6 kpc is 
tfiow — As/u < 3 x 10 7 yrs. However, the radiative 
cooling time in the ambient gas near the radio lobes 
(p « 2.5 x 10~ 26 g cm" 3 , T « 2.5 x 10 7 K) is much longer, 
tcooi « 5m p kT/A(T)pp w 2.7 x 10 9 yrs, where A(T) is 
the usual coefficient for optically thin thermal emission. 
Furthermore, the calculation we consider is embedded 
in the gravitationally supported hot gas atmosphere in a 
galaxy cluster where local cooling after the loss of energy 
by radiative emission does not occur. Any decrease in the 
local gas temperature due to radiative losses is rapidly 
reversed by a global compression and gas inflow that 
maintains the virial temperature and pressure required 
to support the entire atmosphere. Moreover, no ther- 
mal jump associated with the Virgo radio lobe bound- 
aries has been observed in detailed X-ray observations 
(Million et al. 2010). Non-thermal electrons of energy 
7 = E/m e c 2 w (27rm e c 2 /3PA) 1 / 2 « 6 x 10 3 radiating in 



a 3pG field produce synchrotron emission at A = 90 cm 
as in Figure 1. The synchrotron lifetime for these elec- 
trons, t sync = 7/7 = (3m e c/4cr T )(47r/P 2 7) w 1.1 x 10 8 
yrs, also exceeds tfi ow . However, the local synchrotron 
lifetime is not particularly relevant for our diffusion front 
calculation where new cosmic ray electrons are continu- 
ously provided by diffusion from the adjacent radio lobe. 
But we do assume that the synchrotron radio spectral 
index a, appropriate for a power law distribution of elec- 
tron energies, remains constant across the diffusion front. 
Finally, since the flow time across the diffusion front at 
the lobe boundary is almost certainly less than the dy- 
namical age of the lobe itself, the diffusion front structure 
is approximately in steady state during times ~ tfi ow . 

For simplicity, and consistency with current Virgo 
X-ray observations, we consider the cluster gas to be 
isothermal in the vicinity of the radio lobe-cluster gas 
boundary, with pressure, P = pc 2 ., depending only 
on the isothermal sound speed which replaces equation 
(3). The cosmic ray energy density e c is integrated 
over the relativistic energy or momentum distribution, 
e c oc / EN{E)dE oc / p 4 .f(p)(l + p 2 )-^ 2 dp and may 
refer to electrons and/or protons. 

For steady-state, one-dimensional, plane-parallel, 
isothermal flow near the lobe-gas boundary these equa- 
tions reduce to: 

d( P u) n (5) 



dx 







d_ 

dx 



[pu 2 + P + P c ] = 



d(e c u) du d f de c 
dx c dx dx \ dx 



0. 



Equations (5) and (6) are directly integrable 
pu = A and [p{u 2 + c 2 s ) + P c ] = B 



(6) 



(7) 



(8) 



where u < and A < 0. For a steady diffusion front 
propagating in the positive x direction, in our calcula- 
tion, performed in the rest frame of the front, the gas 
flows in the negative x direction from the cluster gas to- 
ward the lobe - the direction of increasing x is shown 
with the small arrow in Figure 1. Eliminating p between 
equations (8) results in a quadratic for u 



+*4(f- B )+'=-= o 



(9) 



for which the desired solution is 




B 



4d 



1/2- 



(10) 



By differentiating equation (9) 

where b = 2Au + [(e c /3) - B], (11) 
equation (7) can be written as 



du 
dx 



u de c 
3b~dx~ 



de c 4e c w de c d ( de c . 
u-r 1 --^ — 7 s --— k— ^ =0. (12) 
dx 9b dx dx \ dx ' 
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The ^-dependence of k can be subsumed into a new in- 
dependent variable 



/dx 
K 



by multiplying each term of equation (12) by k, 

de c 4e c u de c d 2 e c 



d£ 9b d£ dtl 2 



= 



(13) 



(14) 



which can be solved as a pair of first order equations for 
dependent variables <j> = de c /d^ and e c , 



d(j) , A 4e c \ de c 
— =«0^1-— j — 



(15) 



Once e c (£) is known, follows from (10), and the 

density p(£) = A/u(£). 

Solutions of these equations describe a steady state ID 
diffusion front representing the dynamical encounter be- 
tween cluster gas containing cosmic rays inside the radio 
lobe on the small-£ side and approaching external cluster 
gas on the large-.!; side containing no cosmic rays. The gas 
pressure must increase in the ^-direction by P c = e c /3 
to keep the total pressure Pt ~ P + P c nearly uniform 
in subsonic flow, u « c s . Since no X-ray cavities are 
known to be associated with the Virgo radio lobes, it is 
likely that the cluster gas pressure contributes substan- 
tially within the lobes, i.e. the Virgo lobes are partial X- 
ray cavities, only modestly depleted of cluster gas. The 
plane-parallel equations are valid within a distance from 
the lobe boundary that is small compared to the 23 kpc 
radius of the Virgo lobes. The solution profiles in £ - 
e c (£), p(£), u(£), P(£), etc. - are universal, valid for 
any assumed spatial variation of the cosmic ray diffusion 
coefficient n(x). 

Without loss of generality, we assume that £ = at 
some arbitrary position x = within the Virgo radio lobe 
where gas with temperature T = 2.5 x 10 7 K and density 
po = 2.5 x 1(T 26 gm cm~ 3 (n efi = 0.0129 cm" 3 ) flows 
away from the front with velocity uq < and with a cos- 
mic ray to gas pressure ratio (P c /P)o from which e Cj o can 
be found. Our integration of equations (15) begins inside 
the radio lobe at £ = 0, where flow variables have sub- 
script "0", and proceeds upstream toward the diffusion 
front at larger £. The initial gas pressure Po = Pqc 2 . can- 
not be changed much since the density and temperature 
(sound speed c s — 580 km s _1 ) are chosen to approxi- 
mately match those observed in the Virgo cluster near 
the lobes (Ghizzardi, et al. 2004). The velocity uq, cos- 
mic ray energy density e Cj o and its slope </>o = (cte c /c?£)o 
at £ = can be changed to alter the morphology of the 
front. The numerical integration of equations (15) pro- 
ceeds from £ = toward positive £ until e c becomes very 
small. 

We consider two representative subsonic flow velocities 
m at £ = 0. The first is a characteristic velocity for old 
Virgo lobes, uo « —23 kpc/10 9 yr w —20 km s" 1 and 
the second is somewhat larger, uq — —100 km s _1 . The 
upper three panels of Figure 2 show these two solutions 
as functions of the universal ^-independent variable £. 
Notice that as e c approaches zero at some finite £ = £i 
at the leading edge of the diffusion front, the gas velocity 
is asymptotically finite and negative as cluster gas enters 



the front; the second term in equation (14) also vanishes 
at £i, indicating that diffusion at the leading edge of 
the front is balanced by advection downstream. We find 
£i = 2.65 x 10~ 8 and 3.13 x 10~ 8 s cur 1 for w = -20 
and —100 km s _1 respectively. 

For the limiting case of a static lobe-cluster gas front 
where u — > 0, the total pressure P t = P + P c is every- 
where exactly uniform and equation (14) becomes simply 



de c 
d£ \~dj 



0. 



(16) 



The solution in this case is linear e c (£) = e Ci0 + </>o£ 
with uniform slope de c /dt; — cf>o < 0. The gas pressure 
P(0 =Pt~ e c /3 and density p(£) = [P t - (e c /3)}/c 2 pro- 
files are also linear. Since the flow velocities we consider 
| Wo | are much smaller than the sound speed c s , the solu- 
tions plotted in Figure 2 resemble those of static fronts. 
These idealized static solutions are not solutions of the 
time-dependent diffusion equation and represent a math- 
ematical limit that may not be physically relevant. 

Steady state cosmic ray diffusion fronts occur only in 
a gasdynamical environment. In free isotropic cosmic 
ray diffusion in a stationary gas (u = 0) solutions of 
the diffusion equation de c /dt = V-(«Ve c ) depend only 
on a similarity variable £ = r 2 /t and have no steady 
state solution. However in Figure 2 we see that steady 
solutions are indeed possible in a fluid environment. 

2.1. Radio Synchrotron Emissivity 

Computation of the synchrotron emission at a fixed 
radio frequency in the diffusion front requires knowledge 
of the density of relativistic electrons and the magnetic 
field strength and orientation. Radio polarization obser- 
vations of Virgo at 2.8 cm by Rottmann et al. (1996) 
shown in Figure 3 reveal a remarkably high degree of 
linear polarization (> 70%), particularly near the outer 
edges of the lobes where the field direction is highly or- 
dered. Field vectors plotted in Figure 3 indicate that the 
direction of the magnetic field is tangent to the bound- 
ary of both lobes over most of the lobe perimeter. This is 
consistent with the relative compression in Figure 2 near 
the lobe-cluster gas boundary at £ = £i. In our calcula- 
tion the magnetic field evolves passively with the gas but 
does not have sufficient energy density to influence the 
flow parameters. In a one-dimensional flow with approx- 
imately perpendicular field alignment, magnetic flux is 
conserved if B cx p. 

For a power law distribution of relativistic electron en- 
ergies, N(E)dE = KE-PdE with p > 2 the total energy 
density of cosmic ray electrons is 



J E miTI 



EN{E)dE 



K 



(P-2) 



2-p 
rain 



E 2 m - p x ) (17) 



(p-2) mm 



7 2-p 



and therefore K oc e c E^~ n . As the magnetic field slowly 
varies across the diffusion front, the relativistic adia- 
batic invariance of the electron magnetic moment re- 
quires that /3j_7 oc B 1 ' 2 (Sturrock 1994) so E min cx B 1 ' 2 
and K oc e c Bz(P~ 2 \ For electrons with intermediate en- 
ergies, E m i n « E « E max , the synchrotron emissivity 
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from the electron distribution is 

j v oc KB(v/B)- a (18) 

(Longair, 1994) where we assume p = 2.5 and a = (p — 
l)/2 = 0.75. Combining these dependencies, the radio 
synchrotron emissivity observed at the same frequency 
over the whole profile should vary with e c and p as 



where 



I(2p-l) = i(4a + l) 



(19) 



(20) 



so /3 = 2 when p = 2.5. Owen, Eilek & Kassim (2000) 
suggest a mean index a w 1.0 for Virgo for which /3 = 
2.5. The index a may be even steeper near the lobe 
boundaries or vary with x, but current observations do 
not warrant such refinements. 

What is the condition that the emissivity j v oc e c p^ 
(where /3 ~ 2) has a limb-brightened maximum at some 
£ in the diffusion front? Since e c always tends to zero 
as £ — > £i near the leading edge of the diffusion front, a 
maximum is guaranteed if the emissivity increases with 
£ at £ = 0, i.e. 



}_dju_ = ]_^ +/3 l d P, > Q 

j„ d£ e c rf£ pd£ 



(21) 



at £ = 0. Using A = pu and equation (11), the density 
derivative in the second term can be written 



1 dp 1 du 1 de c 
pd£ ~ ~u 14 ~ 3b~dJ' 

The condition for a maximum in j v is then 



e c + 36 



1 + 0- 



Pr 



pu' - pc% 



> 



(22) 



(23) 



at £ = 0. Since </> = (de c /d^)o is always negative, the 
quantity in brackets must also be negative. Furthermore, 
we expect diffusion fronts to be highly subsonic, u 2 << 
Cg, so the condition for a synchrotron emission maximum 
in the diffusion front is 



Pc 

p 



> 



Pr 



4a + 1 



or 



> 



1 



P + P c 1 + /3 4a 



inside the lobe at £ = 0. 



3 

(24) 

Normalized emissivity pro- 
files ji/(0/i^(0) with — 2 are plotted in the lowest 
panel of Figure 2 for the two chosen values of uo and for 
(P c /P)o = 1 which produces the desired maximum for 
radio limb brightening. 

The spatial variation of e c and other parameters with 
the true spatial coordinate x can be found for any as- 
sumed n(x) by integrating dx/dt; = k to determine 
x = x{£). If k is constant, then £ in the solutions in 
Figure 2 can simply be replaced with x/k. 

2.2. Diffusion Fronts in Physical Space 

From any solution in variable £ like those in Figure 2 
an infinite variety of solutions in physical space x{£) can 
be found by choosing various functions n{x) — dx/d£. 
In particular, very sharply decreasing profiles e c (x) and 
j v (x) can be generated even when e c (£) and j v {C) are 



very slowly varying. By this means, x(£) could be cho- 
sen to match the limb-brightened j v {x) and radio flux 
observed at the edge of the Virgo lobes. 

Since the diffusion coefficient n must be positive, the 
spatial coordinate x must be an increasing function of 
£ with positive slope n = dx/d£ > 0. We consider two 
simple examples for x = x(£) normalized so that x(£ = 
0) = 0: 

Case I: — = (1 - e"^ ) (25) 
xq 



dx 



K= dx = ^0 e - mo = 

Co 



Xo 



and 



Case II : — = tanh 
x 

dx x 2 ( O x o 



(I) 



1-1^- 

x 



(26) 

(27) 
• (28) 



For both cases the physical scale of the front extends from 
x = to x — xo where k becomes zero. The parameter £ 
characterizes the rate of distortion of the ^-profiles into 
rc-profiles with smaller values of £0 concentrating the re- 
profiles closer to x = xo 

Figure 4 shows the variation of k(x) and j v {x)/j v (0) 
for both Case I and II with xo = 1 kpc for three values 
of £0, all for the uo = —20 km s -1 solution plotted in 
Figure 2. Figure 5 shows the spatial variation of the flow 
variables on the radio lobe side of the diffusion front, 
p(x), u(x), B(x), P c (x) and P(x) for £ = 3 x 10" 9 s 
cm -1 and xq = 3.08 x 10 21 cm = 1 kpc. The top panel 
shows a rapid decrease in p{x)/po and B(x)/B behind 
the front near x = xo as gas accelerates away from the 
front into the lobe. The decreasing field strength in this 
accelerating flow is responsible for the narrowness of the 
observed limb brightened region at the lobe boundary. 
As x approaches xo from the lobe side, the cosmic ray 
diffusion coefficient k(x) (really k±(x)) plummets (along 
with P c (x)) toward zero. On the cluster gas side of the 
front, x > xo, where e c is assumed to vanish, equation (9) 
requires that the flow velocity u remains locally uniform, 
u(x > xo) = u(xo) along with the gas density and pres- 
sure. In addition, the magnetic field perpendicular to the 
flow continues to be uniform in the cluster gas beyond 
the radio lobe; bright radio rims occur as the diffusing 
electrons encounter this inflowing field. Evidently the 
tangential B field just ahead of the front is prepared in a 
compression between a subsonically outwardly expand- 
ing lobe and/or radiatively cooling cluster gas inflowing 
toward the lobe. 

2.3. X-ray Features Associated with Radio Lobes 

The pressure of relativistic electrons responsible for the 
radio emission observed in Virgo is unlikely to be the only 
or even the dominant non-thermal pressure component. 
The electron energy distribution EN(E) peaks toward 
lower energies that emit at radio wavelengths too large 
to be currently detected. In studies of radio synchrotron 
and inverse Compton X-ray emission from electrons of 
the same energy inside X-ray cavities, assumed to be in 
pressure balance with local cluster gas, Croston et al. 
(2005) find no evidence for an energetically dominant 



COSMIC RAY DIFFUSION 



■5 



relativistic proton population. Although relativistic pro- 
tons dominate the cosmic ray energy density in the Milky 
Way, they are difficult to detect in extragalactic sources. 

Observational efforts should be made to detect the 
total non-thermal pressure in the Virgo radio lobes by 
identifying a lower gas pressure in the lobe compared to 
that in cluster gas outside the lobes at the same dis- 
tance from the cluster center. Small irregularities in 
the azimuthally-averaged radial X-ray surface brightness 
profile in Virgo have been discussed in detail (Churazov 
et al. 2008), but a similar effort should be made to de- 
tect azimuthal variations in the X-ray surface brightness 
and gas pressure that correlate with the radio lobe struc- 
tures. Since the radio lobe boundaries are well-defined, 
it should be possible to detect or set rather tight limits 
on the contribution of non-thermal components to the 
pressure in the lobe. 

To explore this possibility, we constructed an approx- 
imate azimuthally-averaged radial (bolomctric) X-ray 
surface brightness profile in Virgo (using the Ghizzardi 
et al. 2004 observations) and estimated the brightness 
reduction produced by depicting the gas density by two 
inside symmetric radio lobes of radius 23 kpc with cen- 
ters 13.7 kpc from the cluster center in M87, assumed to 
be in the plane of the sky. We found that the X-ray sur- 
face brightness in mid-lobe regions is depleted by about 
0.55 - 0.70 along M87-centered circles of radius 25 - 30 
kpc, circles that also contain large regions of cluster gas 
outside the lobes that can be used for calibration. Near 
the lobe rims the X-ray surface brightness profile is likely 
to have a slope change, not a sharp discontinuity. The 
gas temperature may also differ across the lobe bound- 
ary since thermal conductivity will also be reduced by 
magnetic fields tangent to the boundary. 

However, no lowering of the gas pressure in the ra- 
dio lobes is apparent in the recent Virgo pressure image 
of Million et al. (2010), although low amplitude ther- 
mal pressure variations may eventually be detected with 
dedicated X-ray observations. If the contribution of cos- 
mic rays to the total pressure is lower than we assume 
here, P C /(P + P c ) w 1/3, the condition (24) for radio 
limb brightening requires a larger value of the spectral 
index a. For example, if P c contributes only about 13% 
of the total pressure in the lobes, we require a > 3 for 
limb-brightening, but the structure of the diffusion front 
remains similar to those in Figures 2, 4 and 5. 

The terminology for feedback-related diffuse radio 
emitting regions in galaxy clusters includes mini-halos, 
lobes and emission from X-ray cavities. We stress here 
that all of these extended radio-emitting regions must 
to some extent also be X-ray cavities, although may not 
have been detected with current X-ray observations par- 
ticularly if P c <C P. Another complication that can arise 
in detecting X-ray cavities is the possibility of inverse 
Compton X-ray emission from upscattered cosmic mi- 
crowave radiation, as recently observed in the large radio 
lobes in the Fornax cluster (Tahiro et al. 2009; Seta et 
al. 2011). Additional inverse Compton X-radiation that 
accompanies the radio emission would make detection of 
thermal cavities in the hot cluster gas more difficult. 

3. FINAL REMARKS 

Our simple calculation shows that cosmic rays in the 
Virgo cluster can be strongly confined within certain re- 



gions of the cluster gas. Confinement is aided by mag- 
netic fields that are tangent to the boundaries of these 
regions where the synchrotron emission decreases with 
decreasing cosmic ray diffusivity. Broad gaussian-shaped 
diffusion fronts associated with free diffusion at constant 
diffusivity k in the cluster gas are not present in Virgo 
and cannot be generally assumed. Nevertheless, free dif- 
fusion of synchrotron emitting electrons may occur in na- 
ture and can be detected by spectral flattening near the 
outer edges of gaussian-broadened radio lobes because 
the diffusion coefficient increases for more energetic rel- 
ativistic electrons. 

The radio lobes in Virgo, and those in other clusters, 
are naturally buoyant because the additional partial pres- 
sure of cosmic rays depresses the thermal gas pressure 
(and gas density) in the lobe to match the cluster gas 
pressure just outside the lobe. As these large buoyant 
lobes attempt to slowly rise in the cluster atmosphere, 
they encounter cluster gas that is slowly inflowing due 
to radiation losses, compressing gas throughout the lobes 
and at the lobe-cluster gas boundaries. This compression 
forces the cluster magnetic field to become nearly perpen- 
dicular to the direction of compression, and tangent to 
the diffusion front interface, sharply limiting cosmic ray 
diffusion in this direction. In this sense the propagation 
of cosmic rays in cluster gas is self-limiting. 

Except for the idealized singular case where the mag- 
netic field is aligned precisely perpendicular to the dif- 
fusion front, as the expanding radio lobe and inflowing 
cluster gas mutually compress at glacial velocities, the 
magnetic field is forced to become tangent to the lobe 
boundary, becoming locally perpendicular to the com- 
pression regardless of its initial morphology. Evidently, 
a similar field-orienting compression occurs throughout 
the Virgo radio lobes, as supported by the polarization 
vectors in Figure 3. Such an alignment cannot simply 
be due to so-called "draping" of magnetic lines of force 
along its boundary since it is most unlikely that large 
scale magnetic field lines exist that connect across the 
lobe boundaries and throughout the interior of the lobes. 
Considering the dynamic activity inside the lobes visible 
in Figure 1, the spatial continuity of polarization vec- 
tors in Figure 3 from one observed beam position to the 
next along the field direction almost certainly does not 
trace a single magnetic line of force. Instead, we envi- 
sion a field that naturally aligns perpendicular to a large 
scale one-dimensional compression much like spaghetti 
arranges itself horizontally when served on a dinner plate. 
Magnetic fields of all spatial scales are aligned in this 
manner, accounting for the remarkable consistency of 
the direction and magnitude of the polarization vec- 
tors observed in Virgo. The insensitivity of compres- 
sive field alignment to the spatial scale or local curva- 
ture of the field is particularly relevant to galaxy clus- 
ters like Virgo where the relativistic electron gyroradius 
r g = ~/mc 2 /eB = 5 x 10 12 (7/10 4 )(B/3/iG)- 1 cm is 10 10 
times smaller than the radius of the lobes. The polariz- 
ing magnetic field experienced by synchrotron-radiating 
electrons, when averaged over kpc-sized radio beams as 
in Figure 3, gives the impression of magnetic field lines 
that are coherent over kpc-scales, i.e. that single field 
lines are "draped" . However, the field seen by the radi- 
ating electrons is likely to have a range of different scales 
(radius of curvature) but the electrons responsible for the 
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observed polarized emission radiate from field segments 
that are locally aligned perpendicular to both the line 
of sight and to the large scale compression The radiat- 
ing segments do not need to be (and are probably not) 
connected along a single large scale line of force. 

The region of compressionally aligned B in our dif- 
fusion fronts could extend significantly beyond the lobe 
boundary where j v — >• in Figure 1. A similar cosmic ray 
confinement by boundary-tangent B is also expected in 
X-ray cavities with corresponding radio continuum limb- 
brightening and high polarization near the radio emission 
boundaries at the cavity walls. Detailed radio observa- 
tions of magnetic field polarization and alignment along 
the walls of known X-ray cavities do not appear to have 
been attempted. 

Our steady state solutions describe the upstream dif- 
fusion of relativistic cosmic rays from a radio lobe slowly 
compressing against cluster gas where the cosmic ray en- 
ergy density is negligible. As indicated by radio polar- 
ization observations of Virgo, the magnetic field in the 
upstream cluster gas is aligned largely tangent to the 
lobe boundary and perpendicular to the gas flow. At 
the leading edge of the advancing cosmic ray diffusion 
front, where non-thermal particles first penetrate into 
the undisturbed cluster gas, the diffusion coefficient de- 
clines sharply with increasing magnetic field strength and 
its unfavorable alignment. Nevertheless, the small non- 
thermal pressure introduced by diffusing cosmic rays near 
the leading edge of the front lowers the local gas pres- 
sure and the reduced gas density p = P/(? s drives a gas 
flow that expands back toward the radio lobe in a di- 
rection opposite to the forward advancing diffusion. Un- 
der certain conditions bright radio synchrotron rims are 
expected along the diffusion front at the lobe bound- 
ary as shown in Figure 1. As the steady diffusion front 
is approached from the radio lobe side, the synchrotron 
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emissivity rises because of the increasing magnetic field 
(B oc p) but ultimately declines with the decreasing 
number density of synchrotron-emitting electrons. The 
sharpness of the radio synchrotron emission peaks and its 
decline to zero flux depends on the rate that the diffusion 
coefficient decreases across the diffusion front. Bright ra- 
dio synchrotron rims are expected when the partial pres- 
sure of cosmic rays inside the radio lobe and the spectral 
index are sufficiently large. 

Our analysis in this paper describes lobe rim profiles 
that are in qualitative agreement with the radio image of 
Virgo. As more data become available at several radio 
frequencies, a precise fit to the radio limb brightening at 
the lobe boundaries in Virgo may allow for the first time 
an estimate of the spatial profiles of the field strength 
B(x) and the corresponding cosmic ray diffusion coeffi- 
cient k(x). In addition, it may be possible to determine 
from X-ray observations whether the total non-thermal 
pressure in the Virgo lobes can be estimated from az- 
imuthal gas density and pressure variations across the 
lobes along cluster-centered circular arcs that intersect 
both the radio lobes and cluster gas outside the lobes. 
Indirect detection or non-detection of the pressure of ad- 
ditional unobserved thermal or relativistic particles in 
radio lobes would have broad astrophysical implications. 
Finally, we note that the presence of a gas density jump 
at the lobe-cluster gas transition in Virgo as in Figure 5 
may help to divert and divide the outflowing gas at the 
top of the eastern mushroom-like flow visible in Figure 
1. 
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Fig. 1. — VLA radio image of Virgo at 90 cm (Owen, Eilek, Kassim 2000) showing radio continuum limb-brightening (beam resolution 
is 0.6 kpc). The small arrow at the western rim of the northern lobe shows the direction of increasing x in the diffusion front. 
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Fig. 2.— Two solutions of equations (15) with u = -20 and -100 km s~\ both with T = 2.5 X 10 7 K, p = 2.5 X 10 
(or n e fi = 0.0129 cm -3 ), (P c /P)o = 1, and <?io = —0.01. The second panel from the top shows solutions for e c (£) (dashed-dot lines), P c (£) 
(short dashed lines), P(£) (long dashed lines) and Pt = P + P c (solid lines). Solutions for uo = —100 extend to slightly larger £ than those 
for uo = —20 km s -1 . Bottom panel is the profile j„(£) of the radio synchrotron emissivity normalized to unity at £ = and with f) = 2. 
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Fig. 3. — Polarized intensity map of Virgo at 2.8 cm with superimposed B vector with lengths proportional to the degree of linear 
polarization. Circle at lower left shows beam size (Rottmann, et al. 1996). 
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Fig. 4. — Normalized radio synchrotron emissivity profiles jv(x) / j v (0) for the no = — 20 km s _1 flow in physical space for Case I and 
II diffusion coefficients n(x) (cm 2 s _1 ) both with xo = I kpc. The three profiles in each panel correspond to go = 1 X 10~ 9 (long dashed 
line), go = 3 X 10~ 9 (solid line), and go = 1 x I0~ 8 (short dashed line), all in units of s cm -1 . 
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Fig. 5. — Diffusion front profiles for no = —20 km s~ 1 flow with £o = 3 X 10~ 9 and xq = 1 kpc. Coordinate transformations x(£) for 
Case II are shown in solid and long dashed curves; adjacent dotted lines show profiles for Case I. Upper panel: Gas density and magnetic 
field profiles (solid line) and flow velocity (short dashed line), all normalized to values at x = 0. Lower panel: gas pressure P(x) (solid line) 
and cosmic ray pressure P c (x) (long dashed line) both in units of I0 -11 dy cm -2 . 



